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Phase Diagram of Bilayer Composite Fermion States 
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We construct a class of composite fermion states for bilayer electron systems in a strong transverse 
magnetic field, and determine quantitatively the phase diagram as a function of the layer separation, 
layer thickness, and electron density. We find, in general, that there are several transitions, and 
that the incompressible phases are separated by compressible ones. The paired states of composite 
fermions, described by Pfaffian wave functions, are also considered. 
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I. INTRODUCTION 

Multicomponent fractional quantum Hall effect is of 
relevance to situations when there exists more than one 
species of electron, distinguished by, for example, their 
spin or the layer index. There has been substantial 
interest in multicomponent states of fractional quan- 
tum Hall effect (FQHE) in recent years both because of 
their experimental realizations, and because of theoreti- 
cal progress in the understanding of their structure. The 
focus in this article will be on multicomponent states in 
a bilayer system. For simplicity, wc will assume that the 
Zeeman energy is sufficiently large that the electron spin 
is completely frozen, which is a good approximation in 
the limit of large magnetic fields. 

Multicomponent FQHE states were first considered 
in 1984, when Halperin |0| generalized Laughlin's wave 
function for the v — X/i FQHE |^ to states containing 
more than one component of electron: 

X m',m",m[{zi}, {Wr}] = W {Zj ^ Zk)"^ 
3<k=l 

N2 

X n {Wr - W^r" HiZj - 



r<s = l 



, Ni N2 

exp[--(^|z,P + ^Kn] 



(1) 



3 = 1 



Here, z = x — iy and w = x — iy denote the positions 
of the two components of electrons, labeled by subscripts 
j, k and r, s. The exponents m" and m' are odd integers, 
to ensure proper antisymmetry, and m is an arbitrary 
integer. iVi and N2 are the numbers of electrons in the 
two layers. 

These wave functions were first applied to mixed spin 
FQHE states, i.e., when electrons with both spins are 
relevant. (Such a possibility arises in GaAs because the 
Zeeman splitting is much smaller than the characteristic 
Coulomb interaction energy for typical parameters.) The 
function Xm',m",m is only the spatial part; the full wave 
function is given by 



A[Xm',m",m[{zi}, {Wr}] U1U2 . . .UAfi ^1^2 • ■ -rfAr^ 



(2) 



where A is the antisymmetrization operator and u and 
d are up and down spinors. When the inter-electron 
interaction is spin- independent, the ground state is an 
eigenstate of spin. Only the wave functions with a well 
defined spin quantum number are valid variational wave 
functions. (In principle, due to the spin-orbit coupling, 
other wave functions may also become relevant. How- 
ever, the spin-orbit coupling is weak in GaAs and almost 
invariably neglected in theoretical considerations.) Some 
examples of states that have a well defined spin quan- 
tum number, i.e., satisfy the Fock condition are the fully 
polarized state Xi,i,i &i v — l/i (it is identical to Laugh- 
lin's fully polarized wave function except that the spin 
is pointing in the plane, giving zero z-component, rather 
than normal to the plane) and the spin-singlet state X3,3,2 
at 1/ = 2/5. 

Later, the multicomponent wave functions were also 
considered for bilayer systems. The wave function in 
Eq. ^ also applies to bilayer systems, but with u and 
d playing the role of psewdospinors, indicating the oc- 
cupation of the left and right layers, respectively. When 
the layer separation d is zero, the problem is equivalent to 
that of spinful electrons in a single layer with the Zeeman 
energy set to zero, because the intcr-elcctron interaction 
is pseudospin independent, given by 



V 



2D 



er 



(3) 



independent of the layer index of the particles. In this 
limit, again only a limited class of the above wave func- 
tions has the correct symmetry. However, for finite d, 
the effective two-dimensional inter-electron interaction is 
explicitly pseudospin dependent: For electrons in the in 
the same layer, it is 



but for electrons in different layers, it is 



Vu{r) = 



i\/r^^T(P 



(4) 



(5) 



Here, r is the projection of the distance in the plane and 
t and I denote the two layers. Thus, for finite d, none of 
the Halperin's wave functions can be ruled out for sym- 
metry reasons alone. The simplest non-trivial state is 
the state X3,3,i, which describes an incompressible state 
at 2/ = 1/2. A FQHE at i/ = 1/2 has indeed been ob- 
served in a certain range of layer separation in bi- 
layer systems, believed to be well described by x^,is [§)• 
(Note that the filling factor is defined to be the total fill- 
ing; e.g., = 1/2 corresponds to a filhng factor of 1/4 
in each individual layer.) Another interesting example is 
the QHE state xi^i.i at = 1. It is a non-trivial state 
for a bilayer system. For a single layer, v = 1 QHE 
occurs even for non- interacting electrons, because a gap 
opens up due to Landau level quantization. For the bi- 
layer system, however, the FQHE dX v = 1 originates 
due to interactions, as can be seen by noting that each 
layer has v = 1/2, for which there is no gap for indepen- 
dent electrons. There is experimental evidence for this 
state also (Distinguishing the bilayer v = \ from 
the single layer = 1 is subtle, due to the unavoidable 
presence of some tunneling in the real experimental sys- 
tem, and a systematic study of the evolution of the state 
as a function of various parameters is necessary for its 
unambiguous identification. Further discussion and the 
relevant literature can be found in Rcf. j7|.) Recently, 
there has been a revival of interest in the v ^ 1 bilayer 
state because of experimental observations |8| that have 
been interpreted in terms of a Josephson-like collective 
tunneling 

The physics of multicomponent FQHE states is ex- 
pected to be intimately connected to and derived from 
the physics of the FQHE in single component systems. 
However, Laughlin's wave function, which forms the ba- 
sis of the states discussed above, describes only a sub- 
set of the observed states a,t v — n/{2pn ± 1) in a sin- 
gle layer system. The wave functions in Eq. |l| therefore 
also describe only a subset of a larger class of bilayer 
states. To take a simple example of states not described 
by Xm',m".m, considcr TO = 0, i.e., two uncoupled layers, 
and assume equal density in the two layers. The above 
wave function only applies to FQHE at 2/to', with m' an 
odd integer, even though in reality FQHE will occur at a 
much larger class of fractions given hy v = 2n/{2pn± 1). 
Further, Halperin's wave functions are not applicable to 
compressible states, which, as we shall see, play an im- 
portant role in the bilayer phase diagram. 

The general theory of the FQHE in a single layer is 
formulated in terms of a particle called the composite 
fermion, which is the bound state of an electron and 2p 
quantum mechanical vortices pO| . The states of compos- 
ite fermions are described by the wave function |ic| ] 

V'n/(2pn-Hl) = VlLL'^u Yii^j " Zkf^ (6) 
]<k 

which is interpreted asi^* = n filled Landau levels of com- 



posite fermions, v* being the filling factor of composite 
fermions. Here is the wave function of n filled Landau 
levels of electrons, and the Jastrow factor Y\j^i;{zj — Zk)^P 
attaches 2p vortices to each electron to convert it into a 
composite fermion [Q. The symbol Vlll denotes the 
projection of the product wave function into the low- 
est electronic Landau level. These wave functions are 
known to be extremely accurate representations of the 
actual wave functions. Laughlin's wave function is ob- 
tained for the special case n = 1, and represents one filled 
composite- fermion Landau level. In the limit v — > l/2p, 
the composite fermion filling factor approaches n — > c», 
i.e., the effective magnetic field experienced by the com- 
posite fermions, B* = p(j)Q/n, vanishes. They form a 
Fermi sea here [ p^|jl^ ] , as confirmed experimentally . 
In this limit, $00 becomes the wave function of a Fermi 
sea at zero magnetic field, and il'i/2p describes a com- 
pressible Fermi sea of composite fermions. It is a re- 
markable feature of the composite fermion theory that it 
yields a description not only of all of the incompressible 
FQHE states, but also of the compressible liquids in the 
lowest Landau level. 

The multicomponent generalization of the composite 
fermion theory for the mixed spin states in a single layer 
has been considered extensively in the past ||l^ |l^ . The 
generalized wave functions are 

V'i7(2;il) = TLLL'f n,.n, Zkf' (7) 

where n — ^ ni and Zj label all particles. Here, 
and n| are the numbers of occupied spin-up and spin- 
down composite-fermion Landau levels. This wave func- 
tion satisfies the Fock conditions for all choices of and 
n|, with each choice giving a FQHE state with a defi- 
nite total spin, namely (A^/2)(nf — n|)/(nf +ni). In the 
limit of vanishing Zeeman energy, the ground state has 
the smallest total spin: the total spin is 5 = when n is 
an even integer, with = ni — n/2] when n is an odd 
integer, we have rif = (n + l)/2 and ni — [n — l)/2, with 
S = N/2n, N being the total number of particles. As 
the Zeeman energy is increased, at certain critical values 
of the Zeeman energy ri| increases by one unit while 
decreases by one unit. At very large Zeeman energies, a 
fully polarized state is obtained. This description is in ex- 
cellent agreement with transport and optical experiments 
]l9[ , and also with exact diagonalization studies on small 
systems . It has also been extended to finite temper- 
atures in a Hamiltonian formulation p7| . It ought to be 
noted that the wave functions "01/(2^^1+1) describe only 
uniform liquid states; another structure, a charge/spin 
density waves of composite fermions, has been discussed 
in the context of certain experimental anomalies [^ . 
Halperin's X2p+i,2p+i,2p state at v — 2/{Ap+ 1) is ob- 
tained as a special case with = 1 and ri| = 1; the 
other composite fermion (CF) states are not expressible 



in the form of a Halpcrin wave function. 

Our objective in this work is to generahze the com- 
posite fermion theory to bilayer systems and determine 
the phase diagram of bilayer composite fermion states. 
The generahzation of the composite fermion theory from 
single to two layers is, to an extent, straightforward, as 
we show below. Much work on possible bilayer states 
has been done in the past [^, but these studies do not 
tell us which states would actually occur in nature, and 
in what parameter range. Our goal is to obtain a quan- 
titative description. We determine the phase diagram 
of composite fermion states as a function of two impor- 
tant parameters of the problem, the layer thickness W 
and the layer spacing d (defined to be the distance be- 
tween the centers of the two layers). Due to its varia- 
tional nature, the study cannot rule out other states, but 
we believe that the states considered here exhaust the 
physically relevant uniform liquid states in the parameter 
ranges considered. An important feature of the resultant 
phase diagram is that the incompressible CF phases are 
often separated by compressible phases; i.e., as a function 
of the inter-layer spacing, the transition typically occurs 
from a compressible state into an incompressible state, 
or vice versa. In general, several phase transitions are 
possible at a given fraction. Many of these phase transi- 
tions occur in presently accessible experimental regimes, 
and ought to be observable. 

There has also been exact diagonalization work on two 
layer systems p2|"p^p| . It is restricted to extremely small 
systems, because of the drastically enlarged Hilbert space 
(compared to the single layer problem), and is therefore 
limited in its applicability. The exact diagonalization 
approach is also suitable only for incompressible states. 
Nonetheless, these studies provide indications regarding 
the range of parameters where some of the states, for ex- 
ample z/ = 1/2 bilayer FQHE state, are strongest. Our 
results below are generally consistent with the earlier re- 
sults. We also compare our results with available exper- 
iments at v = 1/2 and v = 1. 

The plan of the paper is as follows. In the next section, 
we describe the generalization of the composite fermion 
theory to two layer system. Section III discusses the de- 
tails of our computational method. In section IV, we dis- 
play the quantitative phase diagram of bilayer CF states 
in the density-layer separation (p — d) plane for several 
quantum well thicknesses at total fillings i' = 1/2, 1/3, 1, 
and 2/5. In Section V our results are compared with 
exact diagonalization studies and experiment. Section 
VI considers bilayer CF states with Pfaflian structure in 
each layer. The paper is concluded in Section VII. 



II. COMPOSITE FERMION THEORY FOR 
BILAYER SYSTEMS 

In the case of a bilayer system, we already know the 
answer in two limits. When the layers are far apart, the 
bilayer state is given simply by two independent single- 
layer composite fermion states. In the other limit, when 
the layer separation d ^ 0, the bilayer problem is iden- 
tical to the problem of spinful electrons in a single layer 
(i.e., pseudospin acts like the real spin), with the Zee- 
man energy set equal to zero. Fortunately, as discussed 
above, the nature of the CF states is well understood in 
this limit. At even numerator fractions, the ground state 
is a (pseudo)spin singlet, and at odd numerator fractions 
it is partially polarized. For the pseudospin singlet state 
the densities in the two layers are equal; when the polar- 
ization is not zero, the densities depend on the direction 
of the total spin, with equal densities obtained when the 
total pseudospin is pointing parallel to the plane. 

To investigate the intermediate situation, we consider 
the following class of wave functions which we expect to 
be particularly favorable energetically: 

*(17i,F2|m) - YliZj - Wr)"V'l7i[{2fc}]V'l7j{w,}] (8) 

where the fully antisymmetric wave function t/V describes 
a lowest Landau level projected single layer state at fill- 
ing factor 77. The state described by '^(v-^.u^lm) will be 
denoted by {pi,T>2\m). When both ?/'i7i and V'i72 are in- 
compressible, {T'i,V2\m) is also expected to describe an 
incompressible state. However, if one or both of ip-^-^ and 
compressible, so is {vi,T'2\m). Halperin's wave 
functions are obtained as special cases for vi ~ \ jm! and 
V2 = l/m". 

The physical interpretation of {pi,T'2\m) is as follows. 
In a single layer system, the CF state is obtained by tak- 
ing an electron state and then multiplying it by a 
Jastrow factor that attaches 2p quantum vortices to each 
electron. In the present case, we start with an electron 
state $sj{zj}]<l>7i2 [{wr}], i.e., a state that has ni filled 
Landau levels in layer 1 and ri2 in layer 2. Now, be- 
cause of the layer index, we have more flexibility in how 
vortices are attached. First of all, because of antisym- 
metry within each layer, we must attach an even number 
of vortices [2pi and 2p2, giving Vi = ni/(2pini -t- 1) and 
V2 ~ ^2/(2^2^2 + 1)] in the relative coordinates of a pair 
of electrons within each layer. However, the number of 
vortices in the relative coordinates of an mter-layer pair 
of electrons can be a third integer, which does not have 
to be even. (The pseudospin part of the full wave func- 
tion in Eq. || takes care of the antisymmetry with respect 
to the exchange of electrons in different layers.) In other 
words, we can attach to each electron two types of vor- 
tices, one seen by other electrons in the same layer and 
the other type seen by electrons in the other layer. 



In the symmetric gauge, the wave function describes 
a uniform state of electrons in two disks, one in each 
layer. The sizes of the disks are determined by the largest 
powers of Zj and w^, which also give the number of single 
particle orbitals in the lowest Landau level in the interior 
of the disk (apart from an unimportant correction on 
the order of unity). The largest power of Zj and Wr are 
Ni/vi + mN2 and Nijvi + miVi, respectively. Here, 
A'^i and Ni are the numbers of electrons in the first and 
second layers, respectively. In order to ensure that the 
electrons in the two layers occupy the same area, A'^i and 
iV2 must be related by 



mTVi 



(9) 



The filling factors of the individual layers are given by 



-1 _ TT-l 



(10) 



(11) 



and the total filling factor is v — i^i + i'2- 

We specialize in the rest of the article to the situation 
when the individual densities in the two layers are equal, 
i.e., Ni = N2. From the preceding considerations, this 
implies Vi —V2 ^ V, vi = 1^2 = v 12, and the total filling 
factor is given by 



2v 



V = Vl + V2 = 



1 



(12) 



Since we are interested in enumerating the states at a 
given total filling factor we write 



(13) 



2 — mv 

Thus, for a given the possible states are 



2 — mv 2 — mv 



(14) 



A. Primary bilayer fractions 

For a single layer, the primary FQHE states t/V occur 
at fractions 



(15) 



2pn+ 1 



We note here that negative values of n are also allowed 
|ll| . These produce the primary bilayer states 



n n 



2pn + 1 2pn + 1 



(16) 



at filling factors 



2n 



V = 



(2p + m)n + 1 



(17) 



The numerator may be either even (for \n\ ^ 1), or one 
(which is possible when |7T| = 1). We will only consider 
the primary bilayer states below. 

States for which the numerator of v is odd and greater 
than one make use of single layer states at fillings V which 
are not primary FQHE states. Consider for example u = 
3/7, for which the possible bilayer states are: 



3 3 
14' 14 







3 3 

TT' TT 



3 3 

8' 8 



(18) 



These states are not likely to occur because the con- 
stituent, single layer, states, do not belong to a primary 
sequence. 

Not all primary bilayer states are physical, though. In 
fact, some of them can be ruled out by noting that the 
inter-layer correlations cannot be stronger than the intra- 
layer correlations for any realistic situation. Which states 
are ruled out is determined by considering the c? — > 
limit, where the state (77, 17|r7i) 



X{{z,-Wrr\{{z,^ZufPi^[{z,}] 



(19) 



must be an eigenstate of the pseudospin. Here, we have 
neglected the lowest Landau level projection, which is not 
relevant to the present discussion because it preserves the 
spin quantum number. As discussed earlier, this state is 
an eigenstate of the pseudospin for m = rric, where 



rric = 2p , n ^ 1 
rric = 2p + 1 , n = I 



(20) 



For n 7^ 1, the state at m = rric is pseudospin unpolar- 
ized ||l^,|l^, whereas for n = 1, the state at m = rric, 
which corresponds to = l/{2p+ 1) is pseudospin polar- 
ized. The latter is just the Halperin Xm,m m states which 
we know to be pseudo-spin eigenstates |2^]. The states 
with m > rric are clearly unphysical, because they corre- 
spond to stronger inter-layer correlations than intra-layer 
correlations. 

Note that a special case is the compressible state at 
ly — 2/{2p + m), corresponding to the limit n = 00. Here, 
for m — rric = 2p, we have v = l/2p, where the pseu- 
dospin unpolarized composite-fermion Fermi sea is ob- 
tained llSlJlq]. (Here, the limit n — > 00 can be taken 
along even integer values of n.) 

The overall picture that we expect on the basis of the 
above considerations is as follows. The value m — {) 
describes far separated, independent layers. The integer 



m is expected to increase in unit steps as the layer spacing 
is reduced, until m = mc is reached. This gives us the 
limit when the two layers are very close. The values of 
m > rric are not physically relevant. 



B. Spherical geometry 

The spherical geometry ||2^,|2^ will be used in all our 
calculations, which is suited for an investigation of the 
bulk properties of the various CF states, due to the lack 
of edges. It considers N electrons on the surface of a 
sphere of radius R, in the presence of a radial magnetic 
field produced by a magnetic monopole at the center. 
A monopole of strength q, an integer or a half integer, 
produces a flux of 2q(f>o {(f>Q = hc/e). The single parti- 
cle eigenstates are the monopole harmonics, Yq^n.m [ p6| , 
given by (with the binomial coefficient (^) is to be set 
equal to zero if either /? > 7 or /? < 0) 



Y 



q.n,rrL 
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i-iY 



n \ ( 2q -\- n 



(21) 



q + n — m 



{v*vy'-'{u*uy , (22) 



where represents the angular coordinates 9 and 
the electron, and 

u = cos(6'/2)exp(-i0/2) 



of 



V = sin(6l/2)exp(i0/2) 



(23) 



(24) 



Here, n is the Landau level (LL) index (n = 0,1,...), 
the orbital angular momcntmn is \q\ + n, and m is the z- 
component of the angular momentum (not to be confused 
with the exponent in Halpcrin's wave function). The 
degeneracy of the lowest LL is 2\q\ + 1, and it increases 
by two in each successive higher LL. 

The composite fermion wave function can be general- 
ized straightforwardly to the spherical geometry. The n 
filled LL state $„ occurs at = {N — n^)/2n. The Jas- 
trow factor Ylj ri^j — Wr) is replaced by Y\j j.{ujVr—VjUr), 
which corresponds to qj = N/2. Noting that the q of the 
product is equal to the sum of the q's of individual fac- 
tors, the state 



(25) 



^ 2pn + 1 2pn + 1 
occurs at Q = mqj + 2pqi + qn, i.e., at 

{2pn + ran + l)iV — {2pn + n^) 



Q 



2n 



(26) 



Indeed, the ratio of the total number of particles to 
the flux, 2N/2Q, gives the correct filling factor v — 
2n / {2pn + mn + 1) in the thermodynamic limit 2N — > 00. 



C. Charge of Elementary Excitations 

We will only consider the ground state wave functions 
in this article, but it is straightforward to calculate the 
charge of the elementary excitation of the incompressible 
bilayer states. Begin with the state in Eq. ^ and add 
one electron to each layer, while holding the flux at the 
value given in Eq. Gu. Now the wave function is given by 



JV+l 



n (^j"^ 



)[PLLL<i>f<^>g']L[VLLL<^l''<i>g']R (27) 



where L and R denote the left and the right layers. The 
q's of the inter and intra-layer Jastrow factors now in- 
crease, so the q for <&, called q', must decrease to 



2n 



(28) 



in order to leave the total Q invariant. At q', the total 
number of states in the lowest n LLs is 2q'n ^ v? = 
N — mn — 2np, which implies that, for -|- 1 electrons, 
there are 2np + mn + 1 electrons in the {n + ly* LL. 
Each one corresponds to an excited composite fermion in 
the full wave function. Since a unit charge was added in 
each layer, the charge corresponding to each composite 
fermion is 



(29) 



2np + mn + 1 



III. COMPUTATIONAL DETAILS 

Despite their simple interpretation, the wave functions 
of composite fermions are too complicated for exact ana- 
lytical treatment. Even for the simplest case, namely the 
Laughlin wave function, it has not been possible to obtain 
any numbers by any analytical method. This is hardly 
surprising, because the composite- fermion wave functions 
describe a strongly correlated state of electrons. (We 
note that the Chern-Simons field theoretical technique 
for dealing with composite fermions |2^,|l^ is not able to 
obtain, from first principles, quantities like the ground 
state energies or the energy gaps. For an approximate 
calculation of such quantities in a Hamiltonian approach, 
we refer the reader to the literature ||29[|.) 

Fortunately, it is possible to obtain exact results for 
the variational wave functions, using a numerical Monte 
Carlo method [^o| . Quantitative numbers are obtained 
without making any approximations, and the results can, 
in principle, be made as accurate as we wish, by running 
the Monte Carlo sufficiently long. In practice, reasonably 
accurate results, correct to 3-4 significant figures, can be 
obtained for up to a total of 2A^ ~ 70 electrons. This 



turns out to be sufficient to extrapolate to the thermo- 
dynamic limit in many cases of interest to obtain results 
that can be compared directly to experiment. 

In a purely two-dimensional system, the interaction is 
given by V^^{rjk) = e\r'^^rk \ ' '^^^re e is the dielectric 
constant of the background material. For this interac- 
tion, the energies scale with e'^/elo, where Ig = y/hc/eB 
is the magnetic length; the only density dependence is 
through the magnetic length. However, realistically, the 
electron wave function has a finite width normal to the 
plane of confinement, which alters the form of the effec- 
tive two-dimensional interaction at short distances. We 
incorporate in our calculations the effects of finite thick- 
ness by working with an effective two-dimensional inter- 
action given by 



e J J [r2 + (zi-Z2)2]i/2 



(30) 



where ^(z) is the transverse wave function. (Here, the 
real quantity z is the distance along the normal direction, 
not to be confused with the complex number Zj denoting 
the position of a particle within the plane.) The func- 
tional form of ^ is determined by self-consistently solving 
the Schrodinger and Poisson equations, taking into ac- 
count the interaction effects through the local density 
approximation including the exchange correlation poten- 
tial l^lj]. The only input parameters are the electron 
density, p, and the shape of the confinement potential, 
which will be assumed to be a square quantum well for 
each layer in this article. The barrier heights of the wells 
are taken to be 276 meV in the calculation, although the 
results are not significantly different from an infinite bar- 
rier. The effect of inter-layer interaction is neglected in 
the LDA evaluation of ^. For further details, we refer 
the reader to Ref. [Q. The self-consistent local-density 
approximation is expected to be accurate on the level 
of 20%. The treatment of the uniform positive back- 
ground charge is somewhat complicated due to the finite 
thickness |]l6| ; fortunately, since we are only interested in 
energy differences, we do not need to consider the neu- 
tralizing background charge explicitly. 

The filling factor is defined to be the ratio of the total 
number of electrons, 2N, to the total flux through the 
surface of the sphere measured in units of (pQ, — 2Q. 
However, for a given wave function, the ratio is 2N de- 
pendent, and has order deviations from its thermo- 
dynamic value, I/. A difficulty arises because the different 
states that correspond to the same thermodynamic v do 
not have the same 2N/N^ for finite systems. In other 
words, they have slightly different densities in finite sys- 
tems. We correct for this by expressing the energies in 
units of e^/ro, where is the interparticle separation, 
and then replacing e^/ro by its thermodynamic value. 
This is equivalent to multiplying the distance r in the 
expression of the interelectron interaction by the factor 



\/ Pn I where pjq \s, the density of the finite system. We 
find that, after correcting for the density in this man- 
ner, the energy differences are only weakly dependent on 
iV and give reliable thermodynamic extrapolations, as 
shown in Fig. (|l|). (In fact, the density correction makes 
only a small correction to the energy differences for the 
large systems considered.) The distance between the par- 
ticles is measured along the chord. 

The energy of the composite-fermion wave functions 
will be evaluated by Monte Carlo method, with the lowest 
LL projection handled by the standard method described 
in the literature |Q . We compute the energy per particle 
to within 0.01% using lO'' Monte Carlo iterations. This 
is done for several system sizes up to 2A^ = 72. 



IV. RESULTS AND DISCUSSION 

We have evaluated the energies of various candidate 
ground states v ~ 1/3, 2/5, 1/2, and 1 as a function 
of iV, and determined the thermodynamic value of the 
energy difference. The energy difference has been cal- 
culated relative to the uncorrelated state 5'(i//2,i//2|o) iii 
each case. We have studied a large class of confinement 
potentials, but will present below results only for three 
quantum well thicknesses as a function of the electron 
density. The distance d separating the layers is mea- 
sured from center to center; its minimum value is there- 
fore equal to the width of each quantum well layer. 

Figs. (H) and (||) show the energy difference IS.E = 
E(T7j!\m) — £'(iy/2,i//2|o) of the various states at several fill- 
ing factors as a function of the layer separation, d/lo, 
with each layer taken to be a strictly two-dimensional 
electron gas. It is clear that all states of the type (F, F|m) 
are in general physically relevant, provided m < rric. For 
ly = 1/3, we have rric = 3; for ~ 2/5, 1/2, 3/7, we 
have rric — 2; and for = I, rric ~ 1. For large rf/Zo, the 
ground state, as expected, is (1^/2, i^/2, |0), and as djlt^ is 
decreased, m increases in steps of unity. For the states 
studied here, the transition occurs from an incompress- 
ible state into a compressible state and vice versa. This, 
however, is not a general feature. For example, the bi- 
layer states at = 4/9 are (2/9,2/9|0), (2/7,2/7|l), and 
(2/5, 2/5|3) which are all incompressible. 

At a fixed filling factor, as d/lo is varied, many tran- 
sitions are predicted to occur, with FQHE disappearing 
for a while and then reappearing. This is reminiscent of 
the case of transitions between CF states with different 
spin-polarizations in a single layer as a function of the 
Zeeman energy. There, however, the transition is always 
believed to be from one incompressible state to another. 

In order to make contact with experiment, we have 
considered the effects of finite width. For each value 
of the finite width, a crossover value of d/lo is deter- 
mined as in Figs. (||) and (||), from which the phase di- 
agram is obtained. The phase diagrams are shown in 



Figs. and (|^). The transitions occur in the 

region c?/Zo ~ 1 — 3, which is an experimentaUy accessible 
parameter range. In recent experiments values as low 
as 1.45 for the ratio d/lo have been obtained by gating 
the sample and varying the density. For zero thickness, 
the phase boundaries in the p — d plane would be vertical. 

It must be stressed that, due to its variational nature, 
our study cannot rule out other phases. For example, the 
bilayer quantum Wigner crystal phase has not been con- 
sidered 1^^. However, we believe that the wave functions 
considered in this work exhaust the likely uniform liquid 
states. 



V. COMPARISON WITH EXPERIMENT AND 
EXACT DIAGONALIZATION STUDIES 

Some of the FQHE states have been investigated in 
the past in exact diagonalization studies p^-|2^,|5t . These 
studies typically deal with a total of 2N — 6 electrons, 
i.e. three electrons per layer; the enlarged Hilbert space 
due to the additional layer index makes it difficult to 
increase the number of particles in the exact diagonal- 
ization study of the bilayer problem. As a result, it is 
not clear to what extent the results are indicative of the 
thermodynamic limit. Also, these studies are not able 
to effectively investigate compressible states, or higher 
order incompressible FQHE states. 

As for the incompressible states, there is general agree- 
ment between our study and the earlier ones, when- 
ever a comparison is possible. For example, the regions 
of stability of (1/5,1/5|0) and (1/3, 1/3|2) in Ref. ||] 
are in rough agreement with ours. At = 1/2, both 
Refs. 1^,^ find that the largest overlap between the 
exact state and the Halperin (1/3, 1/3|1) state occurs at 
around d/l^ ~ 1.5, which lies in the parameter range 
where we find the (1/3,1/3|1) to be stable. At i' = 1, 
our Coulomb result dc/lo ~ 1.2 compares well with 
dc/lo ~ 1.3 obtained in Ref. |^^; they find, however, a 
stronger dependence with the quantum well width than 
we do. (Part of the difference might arise from their 
use of a sine function for the transverse wave functions 
as opposed to our self-consistent LDA.) In some cases, 
however, there are quantitative differences. For example, 
Yoshioka, MacDonald, and Girvin j2^] find (for a pure 
two-dimensional system) that at = 1/3, the (1/5, 1/5 11) 
state is stable in the region 1.5 < d/lg < 4, as op- 
posed to our study that finds the regime of stability to 
be 3.0 < d/lo < 4.0. Or, for the (1/3,1/3|1) state at 
v = 1/2, the overlap does not diminish substantially even 
for d/lo up to 4 |3|,|3- 

The table I compares our theoretical results with ex- 
periments. For 1/ — 1/2, the theoretical value d/lo ~ 3.0 
is in excellent agreement with experiment . (Note that 
the geometry studied in this work is closer to the exper- 
iment of Ref. than of Ref. However, for v = \, 



the theoretical value [d/lo ~ 1.2) is substantially smaller 
than the experimental values {d/lo ~ 1-7 — 2.0). There 
may be several possible origins of this discrepancy. The 
most likely cause is that at smaller values of d/lo, the 
effect of inter-layer tunneling, neglected in theory, is not 
negligible. Indeed, Schliemann et al. |Q find that the 
critical dc/lo increases with the tunneling gap. Also, 
when the two layers are close, the inter-layer interaction 
may also affect the form of the LDA form of ^. 

Another interesting feature of our results is that there 
is a range of parameters, with d/lo < 2, where the v = 
n/{2n + 1) FQHE is stable and at the same time the 
V — 1/2 FQHE also occurs. It may seem that the system 
behaves as a single layer except that it also shows FQHE 
&t V — 1/2. This has been seen experimentally 
However, the nature of the v = n/{2n + 1) FQHE is 
quite different from that in a single layer, which ought to 
be verifiable from the detailed properties of this state. 

Another case of interest \s v — 2/3. Here, the rel- 
evant bilayer states are (1/3, 1/3|0), (1/2,1/2|1), and 
(1, 1|2). The last one needs some explanation. Naively, 
it might seem unphysical, because it appears to have 
stronger interlayer correlations than intralayer correla- 
tions. However, we note that "01 here is not $1, but rather 
Wj^ki^j — ZkY'^^i. In other words, the 1 in (1, 1|2) refers 
to the state n/{2n + 1) with n = —1; it is the v* = —1 
state of composite fermions carrying two flux quanta. 
This is the valid pseudospin-singlet state in the d/lo = 
limit We have not studied 1/ = 2/3 in this work be- 
cause (1, 1|2) involves reverse flux attachment, for which 
our method for treating the lowest Landau level projec- 
tion does not work, for technical reasons. However, it 
appears that while the bilayer state at z/ = 2/3 is incom- 
pressible for the two extreme limits of far separated and 
very close layers, a compressible state is to be expected 
in an intermediate range of layer separation. In the ex- 
perimental studies of Suen et al. ^ ] and Manoharan et 
al. , there is evidence for a direct transition between 
a two component incompressible 2/3 state to a single 
component, incompressible 2/3 state as a function of the 
interlayer tunneling (which determines the gap between 
the symmetric and antisymmetric bands). Our theory, 
on the other hand, neglects any interlayer tunneling, and 
deals only with various kinds of two-component states. It 
is not possible to draw any definitive conclusions in this 
regard from exact diagonalization studies Q . 



VI. PAIRED COMPOSITE-FERMION STATES 

As mentioned earlier the variational nature of our 
study does not exclude other candidate states. To this 
end we consider another class of states that involve intra- 
layer pairing of composite fermions. Such a pairing is 
described by a Pfafhan variational wave function [p8[ , 
studied by several groups in the context of single layer 



FQHE at = 1/2 and = 5/2 [|9|jl6|. The Pfaffian 
wave function at 77 = 1 /2p is written as 



i<j 



(31) 



where P/[M] is the PfafRan of the N x N antisymmetric 
matrix M with components Mij — {zi — Zj)~^, defined 
as 



Pf[M] cx ^[Mi^2M3,4...Mjv-i,jv] 



(32) 



where A is the antisymmetrization operator. Pf[M] 
is a real space BCS wave function and so ip^^ can be 
viewed as a p-wave paired quantum Hall state of com- 
posite fermions. We generalize the single layer paired 
state to the case of a bilayer system by substituting ^p^-^ 
into Eq. |: 



Pf 

(l/2pa/2p|m) 



r,3 



(33) 



These states will be denoted by (l/2p, l/2p|m)''^. It is 
well known that at = l/2p in the single layer case 
the composite-fermion Fermi sea has lower energy than 
the paired composite fermion state in the lowest Landau 
level. |l^;^ Here we check whether or not inter-layer cor- 
relations stabilize the intra-layer paired states and thus 
energetically favor the paired bilayer state. 

We have calculated the thermodynamic energies of 
the possible paired bilayer states at total fillings ly — 
1/2, 1/3, 1, and 2/5 as a function of the inter-layer width, 
d/lo- The crosses in Figs. (^) and (^) show the energy 
of the paired bilayer states measured relative to the Un- 



as a function of 
only those paired 



coupled m — bilayer states of Eq, 
d/lo. [We include in Figs. (||) and 
states which are competitive with the CF bilayer states.] 
On the basis of these studies, we conclude that the paired 
CF states of Eq. ^ are not the ground states at any layer 
separation for any of the filling factors considered. 



VII. CONCLUSION 

We have generalized the composite fermion theory to 
bilayer systems in the absence of interlayer tunneling. 
We have studied a class of wave functions at primary bi- 
layer fractions and made detailed theoretical predictions 
regarding the phase diagram of the bilayer states in the 
density-layer separation plane. One feature that comes 
out of our study is that the phase diagram often con- 
tains alternating stripes of compressible and incompress- 
ible phases. Our results are generally consistent with 
previous theoretical and experimental studies whenever 
a comparison could be made. 
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FIG. 1. The energies of the bilayer states (1/3, 1/3|1) and 
(1/2, 1/2|2) at = 1/2 as a function of l/2Af, where 2N is the 
total number of particles. The energies are measured relative 
to the state (1/4, 1/4|0), with A_E = i?(T7,T7|m) — -E'{i/4,i/4|o) • 
The layer thickness and the layer separation are taken to be 
zero in this plot. In this and in subsequent figures the energies 
are given in units of e^/eZo, where Iq = y' hc/eB is the mag- 
netic length and e is the dielectric constant of the background 
material. 



TABLE I. The experimental and theoretical values for the 
critical bilayer separation dc/lo where the transition takes 
place at = 1/2 and v — 1. 



V 


dcjlo (experiment) 


dc/'o (present theory) 


Reference 


1 


2.0 


1.2 


1 


34 
35 
36 




1 


1.7 


1.2 


1 




1 


1.8 


1.2 


1 




1/2 


3.0 


2.9 


1 






tive to the energy of the uncorrelated state (;^/2, i^/2|0) as a 
fiinction of layer separation. Each layer is taken to be strictly 
two-dimensional in this plot. The various states (F, F|m) are 
explained in the text. 
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FIG. 4. The phase diagram of bilayer CF states a,t v — 1/3 FIG. 5. The phase diagram of bilayer CF states at 

in the p — d plane, where p is the total electron density and d in the p ~ d plane, 
is the distance between the two layers (measured from center 
to center). The phase diagram is given for three quantum 
well widths, lOnm, 15nm, and 20nm. The dashed line shows 
the minimum distance between the layers, which is equal to 
the width of the individual single quantum well layer. 




FIG. 6. The phase diagram of bilayer CF states a,t v — 
in the p ~ d plane. 
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FIG. 7. The phase diagram of bilayer CF states &t v = 1 
in the p — d plane. 



